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We start from a six-band model describing the transition-metal t2g orbitals of half-metallic double 
perovskite systems, such as Sr2FeMo06, in which only one of the transition metal ions (Fe) contains 
important intratomic repulsion Ufb- By eliminating the Mo orbitals using a low-energy reduction 
similar to that used in the cuprates, we construct a Hamiltonian which contains only effective t2g 
Fe orbitals. This allows to treat exactly Upe and most of the Fe-Mo hopping. As an application, 
we treat the effective Hamiltonian in the slave-boson mean-field approximation and calculate the 
position of the metal-insulator transition and other quantities as a function of pressure or on-site 
energy difference. 



I. INTRODUCTION 

In the simple perovskites of formula unit AMO3, the 
transition metal atoms M form a simple cubic lattice. 
In the double perovskites A2MM'0q, this lattice is sub- 
divided into two interpenetrating f.c.c. sublattices, in 
such a way that the nearest transition metal M' lies 
in the sublattice occupied by M and vice versa. A is 
an alkaline earth. Recently, the interest on these sys- 
tems has considerably increased after the finding of low- 
field colossal magnetoresistance (CMR) in SraFeMoOg, 
which remains significant at room temperatureJj The sys- 
tem is a half metallic ferromagnet with a Curie tem- 
perature near 450K. The mechanism leading to CMR 
at low fields is believed to be intergrain carrier scat- 
tering between regions of different orientation of mag- 
netization, which is reduced under an applied magnetic 
field. Other double peroyskites systerns which were-stud- 
ied include BaaFeMoOe,^ Sr^eReOg,! SrzFeWOe,^ and 
the alloy SraFeMoi-^jW^^^OeB SraFeWOe is an insulat- 
ing antiferromagnetic system, and the research on this 
alloy monitors the transition from metallic itinerant fer- 
romagnetism in Sr2FeMo06 to localized-electron anti- 
ferromagnetism in Sr2FeW06. Also, substitution of Fe 
by Co or Mn renders Sr2FeMo06 antiferromagnetic and 
insulatingJj Ordinary ab initio calculations fail to explain 
the electronic structure of Sr2feW06, obtaining a metal- 
lic ferromagnetic ground state.1] This result changes if the 
phenomenological method called LDA+U is usedE 

An issue related with the electronic structure of 
Sr2FeMo06 which is not settled yet is the valence of the 
transition metal ions. This is in principle related with 
the metallic or insulating character, since in an ionic pic- 
ture with oxidation states O"'^, Sr+^, Fe+^ and a closed 
shell Mo"*"^, one might expect that the strong Coulomb 
repulsion at Fe sites Upe brings the system close to a 
Mott transition.El While the Mo-3d chemical shift ob- 
served in optical experiments is practically identical to 
that of M0O3 (indicating formal vj-lence Mo"*"^), the for- 
mal valence of Fe seems to be -1-3.0 This apparently con- 
tradictory result is interpreted by the authors as an in- 
dication of covalency. Mossbauer experiments are inter- 



preted also as indicating covalency.E3 Nevertheless, as we 
shall see, even assuming that all O ions are 0~^, there is 
a certain degree of covalency in the insulating state. Re- 
sults from neutron diffraction experiments obtain majg^ 
netic moments ^pe = 4.1±0.1/iB and ^]\io — OiO.l/is.tiJ 
These values are consistent with an insulating state, sug- 
gesting that Sr2FeMo06 is near a metal-insulator transi- 
tion. 

The metal-insulator transition has been studied theo- 
retically applying the slave-boson mean-field approxima- 
tion (SBMFA) to a six-band model containing the rele- 
vant Mo and Fe t2g orbitals, neglecting the hopping ba- 
tween Mo ions and taking a simplified density of statesH. 
The SBMFA permits to treat Upe — TaJ^ in a way equiv- 
alent to the Gutzwiller approximation.113 

The aim of the present work is to derive an effective 
3-band Hamiltonian, eliminating the Mo sites by a suit- 
able low-energy reduction procedure. This procedure has 
been successfully used in the cuarates. After the origi- 
nal proposal of Zhang and Ricellj that in spite of Cu-0 
covalency, the low energy physics of the superconduct- 
ing cuprates can be described by a one-band model, sev- 
eral systematic studies have derived the different terms 
of this model land used it successfully to calculate several 
properties OlJ The advantages of the effective model is 
that it has a smaller Hilbert space in numerical diagonal- 
izations of finite systems and that the largest interaction 
in the problem Upe (or the on-site repulsion at Cu in 
the case of cuprates) is treated exactly inside an effective 
cell. As a consequence, one expects that approximate 
treatments give better results when applied to the effec- 
tive Hamiltonian rather than the original one. This is 
the case of the SBMFA applied to calculate the metal- 
insulator transition in the cuprates: the results are signif- 
icantlyJjnproved when applied to the effective one-band 
model.Ej 

Our resulting low-energy Hamiltonian contains only ef- 
fective Fe sites with a reduced effective Coulomb repul- 
sion ([/ < 2.8 eV for Sr2FeMo06). It can be described as 
a one-band model containing three pseudospin or "color" 
components describing the xy, yz, and zx t2g orbitals of 
the 3d or 4d shell. For the simpler two-color version 
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FIG. 1. Scheme of the effective hopping between nearest 
dxy orbitals through intermediate O 2p orbitals. 

This is already an effective Hamiltonian in which the 
O atoms, which he in between any two transition metal 



using perturbation theoryH'Eil if late 3d transition metal 
atoms (like Ni or Cu) are not involved. Due to the sym- 
metry of the intermediate O orbitals, one realizes that 
hopping Fe-Mo tpM is two dimensional: it is nonzero 
only between a orbitals lying on the plane a (see Fig- 
ure 1). Thus, it conserves-jColor. We take Upe = 7 eV 
from spectroscopic sudies.L3 Here we derive the other pa- 
rameters of the model for Sr2FeMo06 by fitting_a& ini- 
tio results for the t2g bands obtained previouslyja to the 
corresponding result of Eq. (|l]) in the Hartree-Fock ap- 
proximation. In this approximation U pe is taken as zero 
and epe is replaced by Sp^ — epe + ^Upenp^ , where 
rip^ /3 is the expectation value of fj^fia- (independent of 
i and a) in the Hartree-Fock approximation. We have ad- 
justed the two eigenvalues for wave vector F = 0, and the 
other two for wave vector L = (tt/o, tt/o, it /a), where a is 
the lattice parameter of the f.c.c. structure. This fitting 
gives the values of: eff , Emo, tpM and tMAiicr, a,^a), 
where 7o- lies in the plane a (see Fig. 2). The other two 
independent values of tMAiio', cr', 7) are smaller. For sim- 
plicity, their values are derived using the geometrical rel^ 
tions that correspond to direct Mo-Mo hybridization,c3 
and taking the following relations for the components 
with different angular momentum projection of the. hop- 
ping integrals: (ddS) = 0, [ddir) / {dda) = -0.54.EZl The 
occupation was obtained from the resulting tight- 
binding dispersion, what allows us to derive & £pe from 
^Fe- We have repeated the procedure for several lat- 
tice parameters. Using the numerical derivative of the ab 
initio energy with respect to volume, we obtain the de- 
pendence of the parameters with pressure, as indicated in 
Table I. The most significant change with applied pres- 
sure is the increase in the magnitude of tpM- In contrast, 
£mo — spe decreases with applied pressure. Both effects 
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favor a metallic state by decreasing the probability of 
finding localized electrons at the Fe sites. Note that while 
epe lies ~ 2.5 eV below emo, close to SMp-i This 

is consistent with other ab initio calculations .EI'Lrc3 As a 
consequence, ~ 0.5. However this value is signifi- 

cantly increased, and the amount of covalency is reduced, 
when Upe is treated in a more realistic approximation. 
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FIG. 2. Calculated t2g energy bands (full line and solid cir- 
cles) and fitting using Eq. in the Hartree-Fock approxima- 
tion (dashed lines). The Fermi energy is at eV. The wave 
vectors shown are: W= (tt/h, 0, 27r/a), L= (tt/o, vr/a, 7r/a), 
r = (0, 0, 0), X= (0, 0, 2Tv/a), and K= (37r/2a, 37r/2a, 0). 



«((??)) 


P(GPa) 


tpAiieY) 


iMAf(eV) 


EMo - ffFe(eV) 


16.51 


-25.13 


0.1517 


0.0799 


2.787 


15.71 


-24.95 


0.1955 


0.0956 


2.713 


15.31 


-16.74 


0.2459 


0.1122 


2.623 


15.11 


-9.14 


0.2781 


0.1222 


2.555 


14.91 


-3.40 


0.2960 


0.1296 


2.512 


14.81 


3.45 


0.3150 


0.1361 


2.442 


14.71 


7.04 


0.3242 


0.1390 


2.390 


14.61 


10.72 


0.3338 


0.1433 


2.268 


14.51 


18.15 


0.3510 


0.1531 


2.338 



Table I - Pressure and parameters of Eq. (|l|) for diffcr- 



ent lattice parameters, obtained fifing ab initio energies 
at r and L, with the results of Eq. (|l|) in the Hartree- 
Fock approximation. 

III. THE EFFECTIVE HAMILTONIAN 

In this section, we derive a Hamiltonian containing 
only effective Fe sites, which describes the low-energy 
physics of Sr2FeMo06 or other double perovskites. The 
procedup©! i*-ibasically the same as the one used in the 
cupratesEjcil and consists of the following steps: i) 
change of basis of Mo m.i orbitals to Wannier functions 
centered at the Fe sites, ii) exact diagonalization of the 
cell Hamiltonian Hi at each Fe site i, retaining the ground 
state for each number of particles. These states are 
mapped into the corresponding states of a Hamiltonian 
that contains only Fe t2g states, iii) calculating the in- 
tercell hopping matrix elements {H — Hi) in the re- 
stricted subspace of the cell ground states, and ivlatt. 
elusion of the excited states by perturbation theoryll3'E£l 
We actually neglect the last step because these correc- 
tions are very small, ensuring the validity of the effective 
Hamiltonian. 



A. Change of basis of Mo orbitals 

The Mo Wannier orbitals aia centered at Fe site i are 
constructed Fourier-transforming the mja in the plane a 
at which the orbital a lies, and then transforming the 
dependence on the two-dimensional wave vector back to 
real space: 

" k,l 

With a (A,,) (2) 

where is the number of atoms lying in the plane a, ka 
are two-dimensional reciprocal lattice vectors parallel to 
the plane cr, and A^- are the vectors lying in the plane a 
which connect an Fe site with all Mo sites in this plane. 
The vector A^ does not belong to the Bravais lattice. 
Performing the two-dimensional integral, one obtains: 

4 (—1)"^ — (—1)"" a 
C(A.) . ^^-(^ip^, with A. = -(n.x + n,y), 

(3) 

where riy are two integers of opposite parity and x, 
y are the unit vectors of the plane a. 

Using the inverse of Eq. (^ and some algebra, the 
Hamiltonian in the new basis takes the form: 
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(4) 



with 



CT CT 

+ ?/_Fe ^ flcrfi,afla'fi.cr' + 

(7<(T' 

-^[i7I/(a,0)/,^,a,. + H.c.]. 



(5) 



fia and are vectors of the f.c.c. Bravais lattice and 
the former hes in the plane a. The change of basis has 
introduced hoppings at arbitrary distances, but as we 
shall show, they decay very rapidly with distance. This 
decay is different from those of similar hopping terms in 
the case of the cuprates due to the different nature of the 
Wannier orbitals. The hopping between Fe and new Mo 
orbitals at a distance /x^ = f ('^xX + Tiyy) (where now n^, 
Uy are two integers of the same parity) is: 
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^2 („2_„2)2_2(n2+„2) + i- 



(6) 



The weakest decay is along the diagonal direction = 
Uy. For = Uy — 0, tpM{<y,^) /tpM — 16/7r^ ~ 1.62. 
As a consequence most of the original hopping terms are 
included in J^i^i' ^^^^ '^i^^ be solved exactly. For the 
hopping between different Mo ai orbitals we have to dis- 
tinguish two cases: 



(7) 



where 



(8) 



and the sum over Ao- is restricted to vectors such that 
Ag. + rj lies in the a' plane. This sum is evaluated in 
Appendix A. The sum over rj in Eq. (0) is restricted 
to Bravais lattice vectors such that v — -q are NN in 
the f.c.c. lattice, because otherwise ijvfAf (c, cr' , i/ — rf) 
vanishes. Also, due to symmetry arguments, if tr 7^ cr' 
then tMM{<^,(^' ,1^) = unless v lies in the plane per- 
pendicular to both cr and cr'. The color-conserving 
hopping tMAii'^,'^,'^) retains the same form as in the 
original Mo basis with a rigid shift in the positions of 
the orbitals (from Mo to Fe positions). Some values 



of these hoppings for the experimental lattice param- 
eter are: i_FAf(cr, 0) =0.51 eV, tFM{(J,Sa) =0.17 eV, 
tMAiio', cr, 7cr) =0.14 eV where connects nearest neigh- 
bor ions in the f.c.c. lattice and lies in the plane a. For 
larger distances, tMM is at most of order 0.01 eV. 



B. Diagonalization of the cell Hamiltonian 

Hi can be diagonalized exactly in each subspace of def- 
inite number of electrons and of states which transform 
under the point group operations as a basis function of 
some irreducible representation. . The ground state of the 
most relevant subspaces is mapped into the correspond- 
ing state of a monatomic f.c.c. lattice with t2g orbitals. 
For example, in the subspace of one electron and wave 
functions transforming like cr [xy, yz, or zx), the ground 
state has the form: 



\tla)^iA,fl + B,al)\0). 



(9) 



and is mapped onto the state cj^. |0) of the effective 
monatomic Hamiltonian H^ff- Similarly, the ground 
state for two particles has the form 



|z2aa') =[A2alal,+B2flfl 

+C2{a\jl + flal,m ^ 



(10) 



and is mapped onto the state c\^c\„i |0). In the present 
work and since we are interested in occupations near one 
electron per lattice site, we disregard the states with 
three electrons in the cell. This does not modify the 
results for the metal-insulator transition presented in the 
next section. The inclusion of three-particle states is 
straightforward but complicates the notation and intro- 
duces several new terms in H^ff- For Sr2FeMo06 at 
the experimentally observed lattice constant we obtain 
El = eFe-0-092eV and E2 = 2eFe + 2.442ey, where Er, 
is the on-site energy in H^f / , and the effective Coulomb 
repulsion is U = E2 + Eo - 2Ei = 2.625ey. This strong 
reduction with respect to Upe = 7eV is due to the fact 
that |i2crcr') is dominated by the last term in Eq. (p^). In 
other words, it is energetically more favorable to occupy 
one Fe orbital and one Wannier Mo orbital at that site, 
due to the large intr atomic repulsion in Fe. 



C. The intersite terms 

The hopping terms H — J^i Hi, calculated in the re- 
stricted basis defined above, lead to effective hopping 
terms in Hef / . In addition to the dependence on the 
lattice vector and orbitals involved, the effective hop- 
ping depends on the occupation of the two sites in- 
volved. We denote by ti {i = 1 to 3), the matrix ele- 
ments which correspond to the following transitions: 1) 
\iO,jla) ^ \iW,jO) 2) \ila',jla) ^ \i0,j2aa') 3) 
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\il(Ti, j2a2(J^) < — > |i2(Ti(T4, jlcr2)- The corresponding 
matrix elements turn out to be: 

ti{a,cr',iy) = 2AiA2tFM{o;cr' + A\tMM{<7,<y' 
t2{a,a\y) = [C2{Al + bI) + AiBi{A2 + B2)]t^.i{<j,<j' 

+Ai{AiA2 + BiC2)t^M{(r, 
t3{a, <j', u) = 2(AiC72 + BiB2)(AiA2 + BiC2)t^i{a, a', ly) 

+ (AiCi + BiC2ft^i{a,a' 

where tpM{o', cr' .v) = Saa'tpKiio', v)- The procedure fol- 
lowed so far leads to the following effective Hamiltonian: 

iaa' u 

+t2{(J,a',l'){PaP^+^l + P^oP^+^2) 

+t3{a,a',l^)PaP^+.2} + i^■C.l (11) 

where Pu is the projector over I number of particles at 
the effective site i. Due to the rapid decay of the hopping 
terms with distance, in what follows we shall retain only 
NN effective hopping. At this distance, three hopping 
amplitudes tj can be distinguished for each sum of oc- 
cupations j: a) hopping which conserves color (cr = a') 
and the nearest neighbor vector 7 lies in the plane cr. We 
call this amplitude tjc\\ , b) same as before with 7 at an 
angle of 7r/4 with the plane a. We call it tjc±, and c) 
hopping that does not conserve color (ijnc)- In this case 
for a hopping a ^ a' , "f must lie in the plane perpendicu- 
lar to both a and cr'. The largest of these nine hoppings 
are the • These are the only ones which survive if the 
electronic stucture is approximated as two-dimensional 
for each cr. For Sr2FeMo06 we obtain at the observed 
lattice parameter: fi^ll = 0.197 eV, t2c\\ ~ 0.179 eV, and 
^3c|| — 0.109 eV. The next two in decreasing order of 
magnitude are tine — 0.075 eV and tic± = 0.047 eV. The 
remaining four hoppings lie below 0.03 eV. 

D. Corrections from excited states 

Heff can be systematically improved, including the ef- 
fects of the states neglected in Subsection B by perturba- 
tion theory. The first correction terms are of second order 
in effective hopping terms. Most of them involve matrix 
elements of order 0.1 eV or smaller and denominators of 
order 3eV or larger and can be neglected. The smallest 
denominator occurs for an intermediate state created by 
a term similar to , in which the two-particle eigenstate 
of Hi lies ~ 0.3eV above the ground state. However for 
a non-interacting system, with one electron per site on 
average, the probability of finding a singly occupied site 
is ~ 0.1 and decreases with increasing U. Then, the ef- 
fect of this correction on the energy per site is less than 
O.Ole^ and we can also neglect it for fillings around one 
electron per site. 



IV. THE SLAVE-BOSON MEAN-FIELD 
TREATMENT 

The formalism used here is essentially a generaliza- 
tion to more than t|aip colors (spin up and down) of the 
original formulationlld. The Hubbard-model with orbital 
degeneracy was studied by HasegawaEj extending a pm- 
vious formalism used in thcp^riodic Anderson modelEj, 
and by Fresard and KotliarEll. In our case, the interac- 
tions include correlated hopping which does not conserve 
color. However, the projectors Pu are easily expressed 
in terms of bosonic operators, and the approximation 
remains suitable for our problem. The basic idea is to 
enlarge the Fock space to include bosonic states which 
correspond to each state in the fermionic description. For 
example, the vacuum state at site i is now represented 
as e| |0), where e| is a bosonic operator corresponding to 
the empty site; similarly s|^c|^ |0) represents the simply 
occupied state with spin cr, and so on. The bosonic op- 
erators for doubly and triply occupied states are denoted 
and t} respectively. In this way the projectors can 
be expressed in terms of boson operators (for example 

Pi2 — J2a<a' '^laa'^icrcr') ^^'^ interactions between 
fermions disappear from the Hamiltonian. To restrict 
the bosonic states to those with a physical meaning, the 
following constraints must be satisfied: 

+ ^I'rSia + Y 4aa'diaa' + t\U = 1 

s\aS^'y + Y 4<T<T'diaa' + tjU = cJ^Ci^, (12) 
a' 

To simplify the notation, we introduce the following op- 
erators: 

where 

R - CI - e^e- - V s- , ~ d} d- ,n\-^l'^ 

L.,. = (1 - sl^s,,, ~ J2 A.^'d^^.A^r')-"^. (14) 

with a' ^ a ^ a" . Inside the expression (IH) and usmg 
the first constraint ( p^ ) these operators are strictly equal 
to 1. They are introduced to reproduce the correct non- 
interacti|Hg-|r|esult when the saddle-point approximation 
is made.@l2M~E2l The Hamiltonian takes the form 

i(7 iaa' 

+ E {^o^^o'Ma^G' ,R,,)XIX,,, (15) 

{ij)(T,cr' 

E t2 (a, o' , R,, ) [XIY,.,,, + Yl^^ X,„, ) 

CTl 
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+ Ua, a', R,,)Yl^ B.C.} (16) 

In the saddle-point approximation for the uniform, color- 
independent solution, all bosonic operators are con- 
densed, i.e. replaced by numbers, independent of po- 
sition and color (e| = e, = s, d]^^, = d, = t), and 
their values are obtained minimizing the energy of the re- 
sulting non-interacting fermionic problem under the con- 
straints (p^). For one electron per site in the insulating 
state, one has s = 1/3, e = d = t = 0. In general, 
for a multicolor problem, the condensates for only one 
occupation n are different from zero, and near the metal- 
insulator transition, the values for n — 1 and n -I- 1 are 
infinitesimals of the same order, while the other expecta- 
tion values of the condensates are infinitesimals of larger 
order. Thus, using the constraints we can write near the 
metal-insulator transition: 

e2 = 3^2; s2 ^ i - 2rf2; i = 0, (17) 
o 

where d — *■ 0. The stability with respect to d determines 
the position of the metal-insulator transition. Specifi- 
cally, replacing (^ in (^ the energy up to order d'^ 
takes the form 

E{d)^Ei + {3U + ETB)d^ (18) 

where Etb is the energy of a tight-binding Hamiltonian 
in the f.c.c. lattice, in which the three different NN hop- 
ping amplitudes (denoted by the subscripts ^ = c ||, 
c _L or nc are described as before) are weighted averages 
of the ti^: 

4%/3 4 , , 

Ti = hi + -^hi + (19) 

Clearly, the transition is at the point U = —ETB/i, 
and for the case in which the hoppings do not depend on 
the occupanejj of the sites involved, previous results are 
recovered.c3Eil We have calculated Et b using a mesh of 
286 points in 1/48 of the Brillouin zone of the f.c.c. lat- 
tice. The values of ti^ and U were calculated as described 
in the previous section, and we have looked for the value 
of £ Mo — spe (the most uncertain parameter of the origi- 
nal Hamiltonian) at which the metal-insulator transition 
takes place. The result as a function of pressure is repre- 
sented in Fig. 3. In the same figure we show the value of 
^Mo — spe for Sr2FeMo06 derived as explained in Section 
II from ab initio calculations. If for a given pressure this 
ah initio value lies below the metal-insulator boundary 
obtained with the SBMFA, the system is predicted to be 
a metal. In the opposite case, the mobility of the carriers 
is strongly reduced as a consequence of the electronic cor- 
relations and an insulating behavior is expected. We ac- 
tually see that Sr2FeMo06 is very near a metal- insulator 
transition and that negative pressure can drive it insulat- 
ing, mainly due to reduction of the hopping amplitudes 
as the lattice parameter is expanded. This might be the 



main reason of the insulating character of Sr2FeW06. In 
our formalism, in the insulating phase, the occupation at 
Fe is given by A\. This value at the transition is 0.97 
and should decrease further when delocalization effects, 
not adequately taken into account by the SBMFA, are 
included. Then, one expects an oxidation state near +2 
for Fe in the insulating state, but not exactly Fe^'^. Sim- 
ilarly, the magnetic moment at Fe sites is slightly above 
and that of Mo slightly below zero. Covalency witk 
O atoms, modifies these values, particularly the valence.cl 
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FIG. 3. Full line: phase diagram for the metal-insulator 
transition as a function of smo ~ Spe and pressure. Dashed 
line: values of emo — sfe for Sr2FeMo06, obtained using the 
ab initio results of Table I. 



V. SUMMARY AND DISCUSSION 

Starting from a model describing transition-metal 
t2g orbitals in half-metallic double perovskites like 
Sr2FeMo06, and assuming that on-site Coulomb inter- 
actions are important in only one of the transition-metal 
ions (Fe) , we have derived an effective Hamiltonian H^f / , 
which contains only one species of effective atoms in an 
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f.c.c. lattice with moderate on-site interaction. This 
seems to be the maximum possible reduction of the rel- 
evant Hilbert space per unit cell, and we expect that 
it can be useful for future studies in similar systems in 
which correlations are explicitly taken into account. As 
long as the difference in bare on-site energies (in contrast 
to the Hartree-Fock or ab initio ones) is not very small, 
the interactions in H^ff are important and some treat- 
ment which appropriately takes into account many-body 
effects is necessary for its study. Our derivation can be 
generalized to include both spins. However, our formal- 
ism is not appropriate to treat the case of disorder sys- 
tems in which the highly|-Gprrelated transition metal can 
occupy both sublattices.E3 For an occupation near one 
t2g minority electron per unit cell, Heff is very accurate. 
If the occupation is near two electrons per site, our cal- 
culations should be extended to include states with three 
electrons per unit cell, and perturbative corrections due 
to neglected excited states in the two particle sector. 

We have applied H^ff to analyze the metal-insulator 
transition as a function of pressure or Emo ~ ^Fe i 
Sr2FeMo06. While a similar study was previously done,l 
the accuracy of the present results should be higher be- 
cause the highest energy involved in the problem ([/ pe ~ 
7eV) and most of the Fe-Mo hopping, is treated exactly 
in Heff- We obtain that the region of insulating be- 
havior increases with respect to the previous calculation, 
and a metal-insulator transition is predicted for a neg- 
ative pressure near -5 GPa.. Although we do not know 
the parameters of the starting model for Sr2FeW06, wfij 
believe that the insulating character of the compound,cl 
or the alloys Sr2FeWa:Moi_2;06 for a; ~ 1 is related with 
this transition, ew — ^Fe is expected to he higher than 
£mo — £_Fe due to the larger W-0 hopping.!] Actually for 
X — I, the system is antiferromagnetic and our H^ff de- 
rived here does include both spins. However, we believe 
that the difference in energy between ferro- and antifer- 
romagnetic insulating phases is due to double exchange 
interactions which are smaller than the energetic ingredi- 
ents, whose competition determines the metal-insulator 
transition. In simple terms, in the metallic state, fer- 
romagnetism is favored by energies of the order o£jtlae 
effective tj ~ 0.2 eV (similarly to other proposalsuEj) , 
while in the insulating state, antiferromagnetism is fa- 
vored by energies of order i|/C/ ~ 0.01 eV. The situ- 
ation is similar to the metal-insulator transition which 
takes place in i?Mi03 replacing rare earths R or chang- 
ing temperature.^ It seems that the antiferromagnetic 
order does not affect the metal-insulator transition, and 
the boundary between metal and paramagnetic or anti- 
ferromagnetic order is continuous and smooth. In the 
starting Hamiltonian we have assumed that the O atom 
lies in between its two NN transition- metal (TM) ions. 
If this is not the case, as suggested by some structural 
measurementsH an increase in lattice parameter might 
increase the effective hopping between TM atoms (for 
example-, by a reduction of the TM-O-TM angle, as in 
RNiOsE^). This might be the reason Ba2FeMo06 is still 



metallic in spite of the increase in the distance between 
TM atoms. 
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APPENDIX A: EVALUATION OF THE SUMS S 

Here we describe the evaluation of the sums S{ri) = 
X;a„ C^{A^)Ca'iA„'), with -A„=T]anda^ a'. 
For the sake of clarity, we assume a — xy and a' — 
zx. The other cases are easily obtained using symmetry 
arguments. We can then write: 

f] = ^{n^x + Uyy+n^z); A„j = ^(n^x- n^y); 

A^a; = +nx)x + n;,z], (Al) 

where Ux, n'^, ny, and Uz are integers with Ux + Uy+Uz 
even and n'^ + Uy odd (77 is a vector of the f.c.c. lattice, 
but Act and Acr' are not). Using the definition of A^- 
(Eq.(||)), we can write: 



f NxyNyz 



exp[j — (— fc^n^ + kyUy + k'^{n'^ + Ux) + k'^nz]. (A2) 

We can assume that the system is a cube of edge L, 
and then the number of atoms lying in the plane a is 
Ncr = 2(L/a)^. The sums over the two dimensional wave 
vectors k and k' run over a square Brillouin zone with 
sides of length 2\/27r/a rotated it/ 4 with respect to x. 
Depending on the parity of n^, the sum over n'^ runs 
either over all odd values or over all even values, and: 



E 



exp[j-<(fci, - kx)] = - [Sk^.k' - (-1) 



(A3) 



Using this expression to eliminate the sum over kx in Eq. 
( |A2| ), and replacing the remaining sums by integrals in 
reciprocal space, one obtains: 



27r/a 



'^'('7) = — T / dk' cos(k' nx) / dk' coslk'ux) 

47r-* 7o Jo 



27i-/a-fc^ 



X I / dky cos{kyny) 







-(-1)"" / dky cos{kyny)]. 
Jq 



(A4) 
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The integrals are elementary and after some algebra, the 
result becomes: 



'5' (^^inxX + riyy+n^z 



nj.-\-ny-\-nz ,0 ) 



(A5) 



This expression and Eqs (^) and (0) , allow to calculate 
all hopping terms in the new basis. 
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